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ABSTRACT 

We report the detection of sixteen binary systems from the Anglo-Australian 
Planet Search. Solutions to the radial velocity data indicate that the stars 
have companions orbiting with a wide range of masses, eccentricities and pe¬ 
riods. Three of the systems potentially contain brown-dwarf companions while 
another two have eccentricities that place them in the extreme upper tail of 
the eccentricity distribution for binaries with periods less than 1000 d. For 
periods up to 12 years, the distribution of our stellar companion masses is 
fairly flat, mirroring that seen in other radial velocity surveys, and contrasts 
sharply with the current distribution of candidate planetary masses, which 
rises strongly below 10 Mj. When looking at a larger sample of binaries that 
have FGK star primaries as a function of the primary star metallicity, we find 
that the distribution maintains a binary fraction of ~43±4% between -1.0 to 
+0.6 dex in metallicity. This is in stark contrast to the giant exoplanet distri¬ 
bution. This result is in good agreement with binary formation models that 
invoke fragmentation of a collapsing giant molecular cloud, suggesting this is 


2 J.S. Jenkins et al. 

the dominant formation mechanism for close binaries and not fragmentation 
of the primary star’s remnant proto-planetary disk. 

Key words: brown dwarfs - stars: binaries: spectroscopic - stars: fundamental 
parameters - stars: solar-type - catalogs - radial velocities 


1 INTRODUCTION 


Binary systems, in their various guises, yield vital measures for a range of fundamental 
parameters - mass, radius, luminosity - for the component stars. Studies of the binary pop¬ 
ulation distribution, correlation of orbital elements, and the frequencies of the various forms 
of multiplicity, can be used to shed light on star-formation processes and evolutionary rnech- 


of multiplicity is around 45% 

Raeliavan et al. 2 

o 

1-^ 

o 

perhaps higher than 70% among the 

more massive stars fAbt et al. 

L990; 

Mason et al. 

1998 

Preibisch et al. 19991. and somewhat 

lower (~30-40%) for M dwarfs 

(Fischer & Marc^ 

199/ 

> 


Doppler searches for extrasolar planets have refined the art of single-lined spectroscopic 
binary analysis to the point where relative radial velocities (RVs) of the primaries can be 


-.-l 


trograph stablisation methods (e.g. 

Vogt et al. 

2010; 

Jenkins et al. 

2013d; 

Jenkins & Tuorni 

2014; 

Wittenmver et al. 

2014; 

Anglada-Escude et al. 

2014 

). RV measurements that are not 


tied to an absolute zero point can achieve high internal precision by explicitly removing the 
need to quantify such effects as convective blue-shift and stellar gravitational red-shift; a 
clear account of these and other effects is given in 


Nidever et ah 


(2002). The target stars 


for such planetary searches are generally solar-type and are selected to be chromospheri- 
cally ‘quiet’ in order to minimise the potential for ‘noise’ in any velocity measurement due 


to starspot activity (see 


Jenkins et al. 


2006). They are also selected to have no resolvable 


companions to avoid flux contamination. 

The Anglo-Australian Planet Search (AAPS) selection criteria for its initial group of 


stars are discussed in 


Jones et ah 


( 2002al ). The sample considered here comprises 178 F, G 
and K dwarf stars with declinations south of ~-20° and is complete to F<7.5. There is 
a requirement that the activity seen in CanH&K absorption li nes has an index (measured 
by logi4 K - hereafter referred to as for details see 


Jenkins et ah 


2006 


2008 


2 0111) of 
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less than -4.5, and for there to be no known companions within 2 arcsec. The spectroscopic 
binaries presented in this paper are drawn from this sample. 


2 OBSERVATIONS & DATA REDUCTION 
2.1 The Primaries: Stellar Characteristics 


A summary of the characteristics and masses for the primaries are given in Table CD Spectral 
types, B- V colours, magnitudes, and parallaxes for all the stars are taken from the Simbad 
and HIPPARCOS databases. Metallicities for the stars in this sampl e ar e drawn from two 

02006 ) and photometric 


sources; spectroscopic metallicities were extracted from 


Bond e 


values were taken from the catalogue ofl Casagrande et al. 


al. 


(2 011(1 . For five of the binary stars 


we report there are no Bond et ah spectroscopic values, however, these have Casagrande 
et ah metallicities. In fact, all but one of the primaries, HD39213, have Casagrande et 
ah measurements, allowing for a uniform and self-consistent set of [Fe/H] estimates to be 
generated for the sample. 


data are used with the activity-age relation given in 


Soderblom et ah 


1119911 1 to 

provide secondary age estimates for the stars. Note that since our stars were preselected to 
have f?hi<; values below -4.5, the Soderblom et ah relations are q uantitatively the same in this 
regime to the updated work of Mama i ek fe H i llenbrand (12008! ) due to the sparse activity- 


age data for older c 


warf stars. T h e act iv ity values hav e been draw n from the studies o f 


Henry et al. 

(1996 

), 

Tinnev et al. 

(2002 

), 


Jenkins et ah 


(2006). and 


Jenkins et ah 


(20UJ), 


yet even considering these four works, there are still three stars with unknown activities, 
highlighting the lack of chromospheric activity studies in the southern hemisphere compared 
to the north. The primary stellar ages, along with the stellar masses, are determined through 
interpolation of the Yonsei-Yale isochrones f|Yi et al.l 20011) an d uncertainties can be found 
to reach 100% for the ages of these types of old and Sun-like dwarf stars. Given the ~0.06- 
0.10 dex uncertainty in metallicities, natural variations in stellar activity (for example solar 
R' nK activity variation between -4.75 and -5.10 translates to an age variation from 2.2 


and 8.0 Gyr 


Henry et al.l ). uncertainty in the precise form of the age-activity relationship, 


along with the possibility of flux contamination from the secondary, a number of isochronal 
mass/age/metallicity points can equally account for a star’s colour and magnitude. In fact, 
an offset is found between ages derived from the activity indices and those measured from 
isochrone fitting, whereby the activity derived ages are generally significantly younger than 
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Table 1. The Primaries: Stellar Characteristics 



Star 

B-V 

V 

Parallax 

(mas) 

M v 

Spectral 

Type 

[Fe/H] 

Casa 

[Fe/H] 

Bond 

l°g^HK 

Mass 

(M 0 ) 

1 

HD18907 

0.79 

5.9 

31.1 

3.4 

K2V 

-0.46 

-0.50i0.07 

-5.11 

1.05 ± 0.15 

2 

HD25874 

0.67 

6.7 

38.6 

4.6 

G2V 

-0.02 

- 

-4.95 

1.00 ± 0.05 

3 

HD26491 

0.64 

6.4 

42.3 

4.5 

G1V 

-0.11 

-0.08±0.07 

-4.95 

0.97 ± 0.05 

4 

HD39213 

0.81 

9.0 

16.3 

5.1 

K0V 

- 

0.20±0.07 

-5.10 

0.93 ± 0.05 

5 

HD42024 

0.55 

7.2 

18.2 

3.5 

F7V 

0.19 

- 

- 

1.30 ± 0.05 

6 

HD64184 

0.68 

7.5 

30.0 

4.9 

G3V 

-0.18 

-0.23±0.07 

-4.88 

0.93 ± 0.05 

7 

HD120690 

0.70 

6.4 

51.4 

5.0 

G5+V 

-0.08 

-0.10±0.06 

-4.78 

0.98 ± 0.05 

8 

HD121384 

0.78 

6.0 

25.8 

3.1 

G8V 

-0.39 

-0.40±0.07 

-5.22 

0.98 ± 0.10 

9 

HD131923 

0.71 

6.3 

41.9 

4.4 

G4V 

0.06 

-0.05±0.08 

-4.90 

1.05 ± 0.05 

10 

HD145825 

0.65 

6.6 

46.4 

4.9 

G3V 

0.12 

-0.04±0.07 

-4.74 

1.03 ± 0.05 

11 

HD150248 

0.65 

7.0 

37.5 

4.9 

G3V 

-0.13 

-0.11±0.07 

-4.88 

0.93 ± 0.05 

12 

HD156274B 

0.76 

5.5 

113.6 

5.8 

G8V 

-0.40 

- 

-4.95 

0.83 ± 0.06 

13 

HD158783 

0.67 

7.1 

23.7 

4.0 

G4V 

0.05 

-0.05±0.07 

-4.91 

1.04 ± 0.05 

14 

HD162255 

0.66 

7.2 

24.9 

4.2 

G3V 

0.17 

-0.01±0.08 

- 

1.12 ± 0.08 

15 

HD169586 

0.55 

6.8 

21.4 

3.5 

GOV 

0.32 

- 

-4.92 

1.25 ± 0.05 

16 

HD175345 

0.57 

7.4 

21.3 

4.0 

GOV 

-0.16 

- 

- 

1.05 ± 0.05 


those measured from fitting the isochrones. This result highlights that more work is needed to 
make ages derived from stellar activity relations, or gyrochronology, and those derived from 
evolutionary models, more consistent for old held stars. A consideration of these uncertainties 
enables us to determine a consistent mass range for each star. 


2.2 Determination of Radial Velocities 

2.2.1 UCLES Data 

Observations were made at the 3.9m Anglo-Australian Telescope using the University College 
London Echclle Spectrograph (LTCLES), operated in its 31 lines/mm mode. High precision 
Doppler measurements are made possible by the use of an iodine absorption cell that permits 
detailed modeling of the spectrograph Point Spread Function (PSF). The reader is referred 
to Butler et al. (1996; 2001) for a detailed description, however, the procedure is outlined 
below. 

Multi-epoch spectra at a resolution of R ~ 45,000 are obtained for each star with the 
I 2 cell mounted behind the UCLES, imprinting the stellar spectra with thousands of iodine 
absorption lines in the 5000-6200 A region. Each spectrum can be synthesised from a product 
of a Doppler-shifted copy of a pure stellar spectrum for the star in question (a higher 
resolution stellar template spectrum from which the spectrograph PSF has been removed) 
and an iodine absorption spectrum, all of which is convolved with the spectrograph’s PSF at 
the time of the observation. A least squares fitting process matches this synthetic spectrum 
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with the observed spectrum, and determines up to 14 free parameters (one being the Doppler 
shift, one the wavelength dispersion, and the remainder associated with the detailed PSF 
profile). This fitting process is carried out on 2 A chunks of the spectrum between 5000- 
6200 A and the resulting velocities are weighted by the gradient (<9F/ d\) of the spectral 
profile for each chunk. The mean of these weight ed ve locit ies, c orrected for the Earth’s 


motion relative to the Solar System barycentre 


(McCarthy 


19951) . represents the RV for 


that observation. A barycentric correction is also applied to the Julian dates. The internal 
uncertainty is obtained from the standard deviation of the velocities. This technique has 
demonstrated consistently that 3 ms -1 precision is achievable down to the V—7 .5 magnitude 
limit of the survey for suitably inactive stars over the long term ([Jones et ah 2Q02b). The 
barycentric Julian dates and RV data of our sources are given in Tablc[3] 


2.2.2 HARPS Data 


In order to supplement the velocities measured using the UCLES spectrograph, we performed 
a search of the ESO Archive Facility to determine if any of these stars had high precision 
ESO-HARPS measurements that could be used to increase the phase coverage of our orbits. 
The search revealed that six of targets had been observed with HARPS multiple times, such 
that the inclusion of the velocities yielded much better constraints on the binary orbits. The 
HARPS RV data is also shown in Table[3] 

At this point it is worth briefly discussing the HARPS strategy for measuring precision 
RVs from high resolution and high S/N echellc spectra. The HARPS data is automatically 
processed by the HARPS-DRS version 3.5, with the reduced and analysis quality data on 
the Advanced Data Products page of the ESO Archive website. The act ual reduction and 


analysis method itself is based in general on the procedure explained in 


Bararmgetah 


(19961). 


Unlike the method employed by the AAPS using UCLES, no iodine cell is used by HARPS, 
but instead, precision RVs are measured by maintaining the highest stability possible over 
the long term, but placing the spectrograph in a vaccuum tank to maintain the pressure and 
temperature as stable as possible, and feeding the light to the spectrograph using optical 
fibres. 

The actual RV measurements are not performed in chunks like they are using UCLES, but 


each entire echelle order is used to measure the RV. A weighted binary mask is cons 


that synthetically mimics the position of an absorption line in the star flPepe et al. 


ructed 


2002 ), 
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and a weighted cross-correlation between the stellar spectrum and the binary mask gives the 


RV. The mean of the RV for each order gives rise to the final absolute RV measuremen 


the star, with the uncertainty measured following the procedure in lBouchy et al. ( 2001 ). The 


from 


stability of the spectrograph is maintained by using a calibration Thorium-Argon lamp, that 
is also simultaneously fed to the spectrograph using another fibre, allowin g any drifts in the 
wavelength solution to be measured at the 0.1 ms^ 1 level fjLovis fe Pepe 2007 ). The drift is 
then removed from the measured RV to get the most precise value, however a stab ility of 


less 


2010 ). 


han 1 ms 1 has been found for HARPS data over the long term (e.g. 


Lo Curto et al. 


2.3 Orbital Parameters 

For the analysis of spectroscopic binaries the task is to provide a set of orbital parameters 
(period, P, eccentricity, e, periastron angle, u, time of periastron passage, T p , and the 
projected semi-major axis of the primary, a p sin i) - and a velocity offset, z a , that optimise 
the fit of the equation 
27 t a„ sin i 


z = 


pVT 


: (cos(u + u) + e cos u) + z a 


( 1 ) 


to the n observations of line-of sight radial velocity, z, at true anomalies v (derived from 
the observed times, t, through iteration of Kepler’s equation). A least-squares minimisation 
procedure, invoking several 1DL routines, is used to fit the equation. 

In the search for an orbital sol ution, periods are initially identified via Lomb-Scargle 


periodogram analysis flLomb 


197' 


3 


Scargle 


1982). Orbital solutions are plotted in Fig.Q] 


and summarized in Tablc[2] (i & ii). Where the RV data have a monotonic variation, or 
only one extremum occurs without any clear inflection in the RV variation to constrain a 
second extremum (HD 18907, HD25874, HD26491, HD 131923, HD156274B), convergence 
is reached for a number of different periods. In these cases, where the period is clearly greater 
than the duration of the observations, and where periodogram analysis is least effective, we 
consider only the minimum orbital period in our solutions. Periodograms are shown (Fig. [2]) 
for all the targets where the RV phase coverage is nearly a cycle or more. Where the sampling 
is sparse, aliasing introduces spectral power over a range of frequencies and is particularly 
marked for HD 64184 which has just eight RV o bserv ations and for HD 121384 due to the 
highly eccentric nature of the orbit (see 


O’Toole et al. 


2009). This is perhaps a reminder that 


data sets should comprise more than a dozen observations for the periodogram technique to 
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Table 2. (i) Orbital Parameters 1: Where the RV data extend over more than one cycle or a clear inflection in the RV variation is 
seen to be able to constrain a second extremum, a single set of orbital parameters emerge and are listed below. The quantity asinz 
represents the semimajor axis of the binary system. 


Stax 

P 

days 

U) 

deg 

e 

Tp 

JD- 

2450000 

K 

m/s 

a 

au 

xi 

rms 

m/s 

HD39213 

1309 ± 159 

265 ± 23 

0.2 ± 0.1 

1636.7 ± 2.8 

1265 ± 273 

2.3 ± 0.3 

1.97 

17.1 

HD42024 

76.26 ± 0.03 

347 ± 87 

0.19± 0.01 

1154.8 ± 0.7 

3228 ± 133 

0.38 ± 0.01 

3.22 

11.8 

HD64184 

17.863 ± 0.002 

184 ± 84 

0.249 ± 0.004 

949.50 ± 0.05 

12910 ± 439 

0.130 ± 0.002 

0.93 

2.39 

HD120690 

3800 ± 18 

339 ± 55 

0.34 ± 0.01 

1140 ± 11 

6325 ± 100 

4.73 ± 0.02 

27.7 

15.2 

HD121384 

178.7 ± 0.1 

182 ± 28 

0.84 ± 0.01 

852 ± 1 

10887 ± 2201 

0.61 ± 0.02 

7.19 

15.1 

HD150248 

3272 ± 29 

356 ± 68 

0.67 ± 0.04 

2365 ± 13 

1995 ± 12 

4.36 ± 0.01 

3.19 

3.55 

HD158783 

4535 ± 225 

170 ± 46 

0.05 ± 0.05 

986 ± 27 

2133 ± 527 

5 ± 1 

3.03 

6.62 

HD162255 

48 ± 1 

51 ± 5 

0.06 ± 0.01 

1017.1 ± 0.1 

16223 ± 3419 

0.27 ± 0.01 

2.49 

7.81 

HD169586 

2935± 119 

59 ± 10 

0.4 ± 0.1 

2653 ± 22 

7857 ± 2708 

4.3 ± 0.2 

11.5 

60.2 

HD175345 

312.4 ± 0.1 

277 ± 14 

0.75 ± 0.05 

1256 ± 11 

16099 ± 3675 

0.92 ± 0.01 

3.44 

13.9 


Table 2. (ii) Orbital Parameters 2: Where the RV data have a monotonic variation, or only one extremum occurs without any 
clear inflection in the RV variation to constrain a second extremum, a number of solutions emerge at different periods. For each 
system we list the orbital parameters for the “best fitting” solution having the shortest period. The associated uncertainties are for 
this fit. 


Star 

p . 

1 min 

days 

u 

deg 

e 

T P 

JD- 

2450000 

K 

m/s 

a 

au 

xl 

rms 

m/s 

HD18907 

13770 ± 3528 

314 ± 34 

0.28 ± 0.07 

865 ± 9 

3112 ± 1903 

13 ± 4 

4.31 

9.94 

HD25874 

71108 ± 2308 

274 ± 10 

0.8 ± 0.1 

4382 ± 2453 

2578 ± 253 

33 ± 2 

4.16 

4.67 

HD26491 

9748 ± 1223 

221 ± 38 

0.57 ± 0.05 

250.9 ± 0.5 

4717 ± 2969 

10 ± 2 

3.16 

9.71 

HD131923 

7495 ± 334 

15 ± 79 

0.72 ± 0.02 

528 ± 30 

5784 ± 821 

8.7 ± 0.4 

21.8 

8.60 

HD145825 

6024 ± 163 

149 ± 53 

0.33 ± 0.02 

1095 ± 13 

669 ± 10 

6.54 ± 0.04 

14.4 

7.68 

HD156274B 

191455 ± 13707 

30 ± 85 

0.0 ± 0.1 

5222 ± 10588 

692 ± 55 

63.5 ± 0.2 

11.4 

5.95 


be properly effective ( jHorne fc Baliunasill986l ). The period inferred from the orbital solution 


in each case is marked with a vertical dashed line. Identification of periods is least effective for 
HD 131923 and HD 156274B, where the phase coverage is less than one cycle. Nevertheless, 
the inflections in the RV variation enable a robust Keplerian period to be found. 


3 DISCUSSION 
3.1 The Binary Stars 


HD 18907 : A high proper-motion star classed as K2 V flGrav et al. 


2006) though consistently 


classified as a G5IV prior to this (eg. Evans et al. I1957T) . its colour and magnitude confirm 
that it is evolving away from the main sequence. Chromospherically quiet -5.11, 


Jenkins et al. 


20061), HD 18907 is most likely an old star with an age in excess of 12 Gyr which 


agrees with isochrone fits based on metallicities derived through spectral and photometric 
analysis - for example [Fe/H]=-0.50 dex from Bond et al. and -0.46 dex from Casagrande 
et al. This age/metallicity scenario translates to a stellar mass estimated to be 1.05 Mg. 
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Figure 1 . Radial velocity curves for our stars, with filled circles representing data from the AAPS and filled squares data 
from HARPS. From top left to bottom right we show the stars in catalogue order, HD18907, HD25874, HD26491, HD39213, 
HD42024, HD64184, HD120690, HD121384, respectively. 
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Figure 1 . (cont.) From top left to bottom right we show the stars in catalogue order, HD131923, HD145825, HD150248, 
HD156274B, HD158783, HD162255, HD169586, and HD175345, respectively. 
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HD39213 HD42024 
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HD64184 
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1.00 1.75 2.50 3.25 4.00 

log 10 (period/days) 


HD 121384 
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HD150248 
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Figure 2. Lomb-Scargle periodograms. These are plots of spectral power against log period and are shown for all the stars 
where the phase coverage of the RV data is nearly a cycle or more. Where the sampling is sparse, aliasing introduces spectral 
power over a range of frequencies and is particularly marked for HD 64184. The period inferred from the orbital solution in 
each case is marked with a vertical dashed line. 


Limited phase coverage of the RV data means that the period isn’t well constrained so 
the orbital solution converges to a number of fits with periods upwards of 10000 days. The 
shortest period orbit translates to a secondary minimum mass (M s sin i) of 0.42 ± 0.10 Mg. 
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Figure 2. (cont) 
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HD 25874- Classified as a G2V by 


Gray et al. 


(120061 ) and listed by (IDommanget fe Nvs 


19 94 ) in the Catalogue of Components of Double and Multiple Stars (CCDM), it is identified 
along with a V— 12 companion (position angle 225°, separation 29 arcsec, 1941) as a common 
proper-motion pair. Jenkins et ah derive an R [ IK index of -4.95, suggesting an age of about 
4Gyr, while Casagrande et ah find a metallicity of [Fe/H] =-0.02 dex, suggesting an age of 
9.2 Gyr, in agreement with the subgiant designation. Given the age-metallicity range, a mass 
of 1.00 ± 0.05 Mq is inferred for the primary. Since the phase coverage from both the AAPS 
and HARPS is limited, a variety of solutions exist, however we find the best solution to the 
data gives rise to a very long period low-mass stellar companion, with a period of nearly 
200 years and a minimum mass of 0.33 ± 0.07 M^. 


HD 26491: This is classif ie d as a G1V by 


classified as a G3 ( Houk fe Cowley 


1975; 


Evans 


Gray et al.l (120061 ). and pre viously consi s 


1964J), and is identified by 


Decin et ah 


ently 


( 2000 ) 


as having a Vega-like IR excess. Bond et ah derive a spectroscopic metallicity of [Fe/H] =- 
0.08 dex, which is larger than the Casagrande et ah value of -0.11 dex. An age of 10.5 Gyr 
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is indicated by the isochrone fit which is somewhat older than the 5 Gyr inferred from an 
-Rjjk index of -4.95 (Henry et al.). Accordingly a mass of 0.97 ±0.05Mq is assigned to this 
high-proper motion star. Given the nature of the RV variation (Fig.[U), orbital solutions with 
periods significantly longer than ~ 7000 days are possible, and in all cases the eccentricity 
is greater than ~0.5. An orbital period of 9747.5 ± 1223.2 days (e = 0.57) translates to a 
secondary minimum mass of 0.50 ± 0.15 Mq. 

HD 39213: The Bond et al. metallicity of this K0 dwarf ( Houk 1982 ) is found to be 
[Fe/H]=0.20±0.07 dex. The measured R' UK index of -5.10 relates to an age of 8 Gyr and by 
fitting a range of isochrones with ages from 7.5 and 10 Gyr respectively, a mass range of 
0.93 ±0.05Mq is inferred. Though the orbital parameters appear tightly defined (helped by 
the RV data having a phase coverage close to one cycle), the residuals are relatively high 
(rms=17.1 ms -1 ) and the xl of 1-97 suggests the fit is significant at under the 5 per cent 
level. A Keplerian period of ~1300days and eccentricity of ~0.2 translate to a companion 
with a minimum mass of 0.07 ±0.01Mq - potentially a brown dwarf if the orbit is being 
seen close to edge on. 

HD 42024- Given a Casagrande et al. metallicity of 0.19 for this star, isoch rones rangin g 
in age from 2.5-4Gyr suggest a mass of 1.30±0.05 Mq for this F7 dwarf (Houk 19781 ). 
though the lack of an R' nK index deprives us of a secondary age indicator. Phase coverage 
is broad and the RV values fold convincingly around a period of 76.26 days (Fig. [T]) with 
well defined orbital parameters to yield a companion minimum mass of 0.066 ± 0.003 Mq, 
another possible brown dwarf companion. Residuals are above average with an rms for the 
fit of 11.8 m s -1 . 

HD 64184- This G3 dwarf f Gray et al. 2006 ) is listed in SIMBAD as a variable ( V— 7.49- 
7.55), though no variability flag is marked in HIPPARCOS. The Bond et al. and Casagrande 
et al. metallicities are in close agreement at -0.23 and -0.18 dex respectively, and an age be¬ 
tween 7-10 Gyr is indicated by the isochrone fits, somewhat older than the ~4 Gyr inferred 
from the 7 ?h K index (-4.88; Henry et ah). The stellar mass is estimated at 0.9±0.1Mq. 
The eight RV measurements fold convincingly around a period of 17.86 days (Fig. [I]) and the 
orbital parameters, which are tightly defined by the broad phase coverage, suggest a com¬ 
panion with an M s sin i of 0.170 ± 0.001 Mq and an orbital separation of 0.130 ± 0.002 AU. 
With a sufficiently high orbital inclination (^87°), the secondary could provide sufficient 
obscuration of the primary’s surface for a variation in V of ±0.03 m agnitudes, given that a 
mass of 0.160 Mq would be no brighter than My ~ 9.6 at 6 Gyr - (IBaraffe et al.l 1998S ). 
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HD 120690: This G5+V, as classified by Gray et ah (200a), has an index of -4.78 
(Henry et al.), equivalent to an age of ~3Gyr. Metallicities are found to be -0.10 and - 
0.08 dex from Bond et al. and Casagrande et al. respectively, suggesting an age from the 
isochrone fits somewhere between 6.5 and 9 Gyr, from which a stellar mass of 0.98 ± 0.05 Mg 
is inferred. While limited phase coverage of RV data means that the orbital period is only 
poorly constrained by periodogram analysis (1900 ± 500 d), a full orbital solution yields 
tightly defined parameters. A Keplerian period of 3799.98 ± 18.12 days and an eccentricity 
of 0.34 indicate a secondary minimum mass of 0.59 ± 0.01 Mq. With the secondary con¬ 
tributing upwards of ~0.4 per cent of the flux at 5500 A, the central region over which radial 
velocities are determined, contamination of the primary spectrum could be a contributory 
factor. It could also affect the star’s colour-magnitude location making it appear redder (and 
apparently more evolved) possibly accounting for an element of the age discrepancy men¬ 
tioned above. These issues are considered further in Sect ion [3~ll We do note that a second 
signal may be present in the data, with a period of 531 days and a semi-amplitude of 32 m 
s -1 , which if it were a genuine Doppler signal, would give rise to a planet with mass around 
lMj. The addition of this signal can serve to decrease the rms by a factor two. 


HD 121384: Classified as G8V fjGrav et al . 20061 ). and listed in iDommanget & Nvsl ( 19941 ) 
along with a common V —13 proper- motion companion (pos ition an gle 45°, separation 


31 arc.• sec. 1941), this star is identified by lOudmaiier et al. (19921 ) and by 


Alima,nn b. Probst 


(119911) as having a Vega-like IR excess. Bond et al. find a [Fe/H] value of -0.40±0.07 dex, 
in excellent agreement with the value found by Casagrande et al. of -0.39 dex, yet in good 
agreement within the uncertainties with Bond et al. The colour and magnitude are well 
fit by 7-10 Gyr isochrones indicating that the primary component is indeed evolving away 
from the main sequence - a view further evidenced by the low level of /?/ HK activity (-5.22, 
Henry et al. 19961 ). Accordingly its stellar mass is estimated at 0.98±0.10Mq. The 179-day 
period is sharply defined by the broad phase coverage and the Keplerian solution, albeit with 
a relatively large rms of 15.1 m s -1 . The data indicate there is a 0.17 ± 0.01 Mq companion 
orbiting with an eccentricity of 0.84. An eccentricity of this magnitude (Fig. HD places this 
binary in the extr eme upper tail of the eccentricity distribution for systems with periods 


< 1000 days ( Dnqnennov fe Mayor 1991), Fig. 6a) and when combined with the window func¬ 


tion, makes the detection of this signal from periodogram analysis alone very difficult (see 
Fig. 12). 
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HD 131923: The colour-magnitude location is well fit by 9-12 Gyr isochrones with metal- 


licities of -0.05 and 0.06 dex for Bond et al. and Casagrande e 


solar values suggest that this high proper-motion G4V star flGray et al. 20.06) is starting to 


al. respectively. These close to 


evolve away 

iom 

Henrv et al. 

1996) 


: rom the main sequence; by contrast, the age inferred from its R ' uK index (-4.90, 


19961) is only ~4 Gyr. Given the monotonic variation in RV measurements, the 


period is poorly defined, though a 7496 day Keplerian solution emerges with an eccentricity 
of 0.72. The rms for the fit is above average (~9m s -1 ) and this is reflected in a xt °f 21.8. 
With an inferred mass for the primary of 1.05 ± 0.05 Mq, the orbital parameters translate to 
a secondary minimum mass of ~0.52 ± 0.06 Mq. The star is identified in the HIPPARCOS 
catalogue as a ‘suspected non-singl e’ object. 


HD 145825: This G3 dwarf (Torres et al 


2006) again appears to have a metallicity con¬ 


sistent with solar ([Fe/H]=-0.04 and 0.12 dex from Bond et al. and Casagrande et al.) and 
isochrone fitting suggests a n age un der 3 Gyr, which is consistent with the relatively high 
level of activity (-4.74, 


Henry et al 


19961 ). Consequently the stellar mass is estimated 


at 1.03±0.05 Mq. Recent RV measurements have improved the phase coverage, leading to 
a more sharply defined period of 6024 ± 163 days. With a minimum mass of 0.06 ± 0.01 Mq, 
this adds to the list of brown dwarf candidate companions in our sample. 


HD 150248: Consistently classified as a G3V (Gray et al. 


2006 


Houk 


1978 


Evans et al. 


19571 ). the sub-solar metallicity of [Fe/H]=-0.11 dex from Bond et al. and -0.13 dex from 
Casagrande et al. suggest an age from isochrone fits of ~ 7 Gyr, which contrasts with 
an activity-inferred age of ~4Gyr (R ' nK activity of -4.88 from Henry et ah). A mass of 
0.93 ± 0.05 Mq is estimated for the star. The orbital parameters are fairly well defined, 
with a period of 3272.4 ± 28.7 days and eccentricity of 0.67, relating to a companion with a 
minimum mass of 0.10±0.02Mj. 

HD 156274B: Listed as a multiple-star system in Dommanget & Nvs ( 19941 ). comprising 


four known components: G1666A, a G8V (Eggl 2013), G1666B, an M0 dwarf (separation 
7.5arcsec, 1880), a V=12.5 companion and a 14.0 companion (respectively 279°, 41.8arcsec, 
1900; 30°, 47arcsec, no year). Our RV measurements indicate that G1666A is itself a spectro¬ 
scopic binary. No spectral metallicity was determined but Casagrande et al. find a metal-poor 
value of [Fe/H]=-0.40 dex. Such a value demands isochrones of 9-10 Gyr, though the star 
would appear to be no older than ~6 Gyr judging from its index (-4.95, Jenkins et al.) 
in which case isochrones with metallicities of 0.1-0.2 dex provide complementary fits. These 
two scenarios translate to a stellar mass estimated at 0.83±0.06Mq. The best fit period is 
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found to be 524 years, implying a secondary minimum mass of 0.20 ±0.03Mj, the fit has an 
rms of 5.95 m s^ 1 , however we note that the uncertainties on the measured quantities from 
the bootstrap are formal to the solution presented, ft is clear that this is the minimum best 
fit to the data and so the true solution could be very different and therefore properties like 
the time of periastron passage are essentially unconstrained. 

HD 158783: This G3/G5 dwarf (Houk & Cowley 1975 1 appears to have a metallicity 
slightly under solar (-0.05 dex from Bond et al. and 0.05 dex from Casagrande et al.) and 
isochrone fits demand an age o f 8-10 Gyr, significantly above the 4 Gyr age inferred from the 
star’s index (-4.91. Hen ry et ah 19961 ). A stellar mass of 1.04 ± 0.05 Mq is inferred. Or¬ 


bital parameters are tightly defined given the limited phase coverage of the RV measurements 
and the Keplerian fit is significant at the 20 per cent level. A period of 4534.78±224.92 days 
and zero eccentricity translate to a secondary minimum mass of 0.20 =ly0.02 Mq. 


HD 162255: Though no B\\ K index is given for this G3 dwarf flHouk fe Smith-Moore 


19881), isochrone fits based on ~-0.01dex (Bond et al.) and 0.17 dex (Casagrande et al.) 


metallicities suggest an age ranging from 5-9 Gyr, translating to a stellar mass of 1.12 ± 0.08 Mq. 
The eleven RV measurements fold convincingly around a 47.95-day period, indicating a com¬ 
panion minimum mass of 0.333 ± 0.001 Mq. 

3) appears metal rich ([Fe/H] = 0.32 dex, Casagrande 


HD 169586: This GOV star 


Houk ( 1982 


et al.) suggesting an age around 2-4 


Gyr, which is in agreement with that derived from the 


i?(j K index (-4.92. iHenrv et al.l ll996h; the stellar mass is estimated at 1.25±0.05Mq. The 


acquisition of two of the most recent RV measurements has defined a relatively sharp ex¬ 
tremum in what was originally a monotonic RV variation so that a 2935-day orbit with an 
eccentricity of 0.35 appears well constrained. Given that the secondary has a minimum mass 
of 0.68 ±0.22 Mq (equivalent to an My = 9.4 companion) contamination of the primary’s 
spectrum has almost certainly taken place, and is possibly a reason why the rms for the fit 
(60.2 m s -1 ) is so large, along with the large Xv- Nevertheless the existing RV measurements 
clearly indicate that the primary has o ne or more companions. 


HD 175345: Listed in 


Dommanget fe Nvsi I1994J as having a V=14.2 proper moti on com- 


panion (B413: 252°, 5.4arcsec, 1927), our measurements indicate that this GO dwarf (Houk & Smith-Moore 


1988) is itself a spectroscopic binary. There is no /?/ HK index for this star and the isochrone 


fits based on its Casagrande et al. metallicity of [Fe/H] =-0.16 dex gives an age range of 
4-9 Gyr along with an inferred mass of 1.05 ± 0.05 Mq. There is a clear fold of the RV data 
around a 312-day period so orbital parameters are tightly constrained, however the residuals 
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Figure 3. (a) The A APS binaries are indicated on a colour-magnitude diagram constructed using HIPPARCOS data for 
solar-neighbour hood stars, (b) Best-fitting isochrones for each star; the numbers assigned to the stars are given in Table[I] 


are relatively large (rms = 14 m s^ 1 ) making the fit statistically poor, even though the Xv 
is fairly low (3.44). With a minimum mass of 0.48±0.08Mq, (equivalent to My = 9.5 or 
brighter) the secondary is again a source of spectral contamination. 

The locations of these stars on a HIPPARCOS-based HR diagram, along with a summary 
of the isochrone fits, are shown in Fig.[3]a&b, respectively. We also show the distribution of 
eccentricity versus period in Fig.[4j We find that the companions to HD121384 and HD175345 
have very high eccentricities for companions with orbital periods below 1000 days, placing 
them in the extreme upper tail of the distribution in this parameter space. We also find a few 
of the longer period companions that have not been fully constrained yet due to the limited 
baseline of the data have moderate to high eccentricity. Although the periods of the orbits 
could be significantly longer, the eccentricities are rather well constrained with the current 
data in hand, assuming a single Keplerian fit. The one obvious exception is HD156274B 
which currently only exhibits a linear trend with velocity over time, and although we fit this 
will a circular model, this could very well be highly eccentric. 


3.2 High Contrast Observations 

A number of the binaries we have discovered in this work have been followed up using adap¬ 
tive optics systems to search for direct confirmation of the companions. Obtaining orbital 
motion from direct images of low-mass companions to bright Sun-like stars, espe cially whe n 
coupled with RV information, can yield dynamical masses for the companions (jLiuj 200711 . 
Combining dynamical masses with photometric colours and spectra can allow evolutionary 
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Figure 4. A plot of eccentricity versus period for the AAPS binaries. The binaries represented by filled circles are from 
Table|2](i) and have well constrained Keplerian periods; open circles are used for binaries with poorly constrained periods - 
given in Table[2](ii)- 


a nd a tmospheric properties to be well constrained and models to be tested (e.g. 


Dupuy et al. 


2009). 


In 


Jenkins et al. 


(120101) we observed two of the host stars we report new binaries for 


in this work, HD25874 and HD145825. In that work we found contrast ratios of greater 
than 11 magnitudes at separations of only 0.5" using the VLT NAOS CONICA instrument 


(IRousset et al. 


2003 ) in Simultaneous Differential Imaging mode. Although a tentative de¬ 


tection of the companion around HD25874 was discussed, further analysis revealed this to 
be a probable artifact of the reduction and analysis procedure, and therefore no companion 
detection was conclusively made for either of these stars with rnid-T dwarf masses of around 
50 Mj or so. Some of the other stars we report companions for are included in our ongoing 
NACO/NICI imaging campaign. 
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HD175345 



JD-2451000 

Figure 5. An example astrometric plot of the positional residuals and associated errors for HD 175345. 


3.3 HIPPARCOS Astrometry 


Out of these sixteen binaries, four (HD 39213, HD 120690, HD 121384, and HD 175345) have 
well constrained periods in the range 0.5-6 years, making them suitable candidates for anal¬ 
ysis of their HIPPARCOS astrometry. In order to determine if there are any astrometric 
signatures that would allow us to place additional constraints on our orbital parameters, 
we extracted the astrometric data from the HIPPARCOS database ( van Leeuwen 2007) and 
derived positional residuals (a cos <5, S). In Fig.[5]we show these with associated errors for 
HD 175345. 


There are no significant variations obvious in the plots for the four candidates, though in 
the case of HD 175345, periodogram analysis reveals significant power at around 320 days, 
(c.f. a 312-day period from RV analysis). Failure, however, of these data to fold around either 
period suggests that the spectral power owes more to sampling of the data (which have a 
strong 300-day element) than to any possible astrometric signature. The significance of the 
astrometric variations is clearly low and we do not attempt further to constrain our orbital 
parameters using these data. 


3.4 Secondary Flux Contamination 

In those cases where the orbits are fully constrained and the secondary minimum mass 
is ~0.5M o (HD 18907, HD 26491, HD 120690, HD 131923, HD 169586 and HD 175345) the 
contribution of the secondary to the overall flux is at least 1 per cent and the signature of this 
contribution is found in the quality of the orbital solutions: the rms scatter of these stars is 
generally much higher than the internal measurement uncertainty. The reason for this is that 
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flux from the secondary, which is contaminating the primary’s spectrum, will be associated 
with a different radial velocity at each of the subsequent observations from that when the 
template spectrum was acquired. The radial velocity fitting process relies on the assumption 
that the primary’s spectrum is modified only by the Doppler shift of the primary and the 
spectrograph PSF variation. Consequently, less than optimal solutions can be expected when 
a faint secondary contribution to the primary’s spectrum is present at variable Doppler shifts. 
In essence, the fitting process matching the observed and synthetic spectra would be expected 
to generate larger measurement errors, and this is the case particularly for HD 169586. 
Fig. [6] shows how the quality of the orbital solutions - measured by the significance of the 
fit - varies with the number of observations (plot a) and the binary minimum-mass ratio 
(plot b ). While no correlation appears to be evident between the reduced Chi-squared and 
companion mass ratio, in general, the statistical significance of the fits appear to indicate 
a marked deterioration as the number of data points increase. This could reflect the fact 
that these are the systems exhibiting the most stellar flux contamination and therefore they 
required more observations in order to better constrain their orbits, or, another possibility 
is that these systems also contain additional companions, either brown dwarfs or planets, 
that are giving rise to mixed signals that are being manifest once enough RV measurements 
have been acquired. 


Given the precision with which radial velocities can be obtained by the AAPS, we could 
speculate that where low Msin i values yield statistically poor fits (HD39213, HD 121384), 
the orbital inclinations are low and that the companions have masses high enough for spec¬ 
tral contamination to be taking place. The large errors in RV measurements for HD 39213 
and HD 156274 lend credence to this scenario, though this isn’t the only explanation. It is 
also possible that these are multiple systems for which a double-star solution simply isn’t 
appropriate. Moreover, enhanced activity in any of these stars (bearing in mind that the 
A' HK index is merely a ‘snapshot’ measurement) would mean that the internal errors would 
not properly reflect the uncertainty in the RV measurements. This is why it is pragmatic to 
use relatively inactive stars in Dopple r sear ches. For the targets considered here, enhanced 

(119961 ) estimate that 90 percent of the time 


activity is unlikely to be the cause: 


Henrv et al. 


a single R ' uK measurement for a solar-type star is sufficient to identify correctly if it has an 
activity greater, or less, than Ah K = -4.75, and all but one of the stars have A^k indices lower 
than this. Where the measurement errors are low but the orbital solutions are statistically 
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Figure 6. Quality of orbital fit (in terms of the statistical significance attributed to the xi value) versus (a) the number of 
observations per star (b) the minimum mass ratio of the system. 





Figure 7. Effect of binarity on the colour and magnitude at 5.5 Gyr for solar metallicity primaries ranging in mass from 
0.88-1.20Mq. Secondary mass values increase in 0.05M P increments from 0.5-1.0M p . 


weak (HD 121384 and HD 162255), the most convincing explanation is that these systems 
comprise more than a single companion. 

The effect of binarity on the colour and magnitude of a primary of any age and metallicity 
can be modelled easily using the Yonsei-Yale isochrone data. Fig.|T| shows the variation in 
colour and magnitude for a 5.5 Gyr, solar metallicity isochrone ranging in mass from 0.88- 
1.20 Mq (solid curve). The effect of adding a secondary is marked by the dotted curves at 
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various mass intervals, with each point representing a secondary companion that increases 
in mass from 0.5-1.0 M p in steps of 0.05 M p . 

For companions below 0.55 M p the effect of binarity on the colour-magnitude location 
of the primary is negligible. The effect of a 0.70 M p companion is to make the primary 
appear redder and brighter by around 0.05 and 0.3 magnitudes respectively, for a solar- 
mass star, with slighter larger and smaller values for less massive and more massive primary 
stars, respectively. Such an effect generally makes the unresolved pair appear older and/or 
more metal rich; however if the primary itself is evolving away from the main sequence, the 
companion will make the pair appear ‘bluer’. These effects will complicate the process of 
mass estimation. For example, the presence of a 0.70M p companion translates effectively to 
a systematic error of ~ 0.04Mq in the mass of the primary. This is of the order of the error 
in mass due to the age-metallicity uncertainties. The orbital solutions indicate binary mass 
ratios generally significantly less than unity (ffD 169586 being the exception) so that the 
effect of binarity on the determination of the primary mass (and by extension the secondary 
mass) is negligible. Clearly, for many of the binaries in our sample, the uncertainty in 
secondary mass is due principally to poorly constrained orbital parameters. 


3.5 Mass Distribution 

The binaries and planetary companions to solar-type stars reported by the AAPS provide 
an opportunity to examine, from a single radial velocity survey, the distribution of Msin i 
values for a mass regime extending from Jovian through brown dwarf to sub solar in value. 

In order to derive a more meaningful distribution of minimum-mass ratios, we need to 
impose a period cnt-off on both the planetary and binary companions so that we count 
only those companions within a certain distance of the primary stars. As a rule, plane¬ 
tary candidates are announced when the phase coverage of the RV data are close to one 
orbital period. For the AAPS, which has been operating since 1998, we can say that the 
inventory of exoplanet candidates orbiting with periods up to 12 years (i.e. out to around 
the orbit of Jupiter in the solar system) is reasonably complete down to the level permit¬ 
ted by a Doppler precision of ~3ms _1 , i.e. complete for Jnpiter-mass objects and above. 
Brown dwarf and low mass stellar companions induce larger reflex velocities making them 
easier to detect over a greater range of distances and periods. In order to compare directly 
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log(Companion Minimum Mass/Mj) 


Figure 8. The raw distribution of companion minimum masses uncovered by the AAPS with periods less than 12 years. 


with 


Du a uenn ov & Mavori (119911 ). six out of the sixteen binaries have periods greater than 


12 years and these need to be excluded from our count. 

Ten remaining binaries out of a sample of 17 8 stars is around half of that expected from 
the period distribution given in Fig. 13 of Ragh a van et ah (2010), when normalised by the 
multiplicity fraction, though there are several reasons for this. First, our sample excludes 
all known short-period spectroscopic binaries with separations less than 2 arcseconds; sec¬ 
ond, all binaries beyond 2 arcseconds detected using Hipparcos data are excluded from our 
sample; third, ‘double-lined’ spectrum binaries are immediately removed as soon as they 
are recognized; fourth, our requirement for the AAPS target stars to have an R ' uK index 
below -4.5 has the effect of filtering out some short-period, chromospherically active bina¬ 
ries - RSCVn/WUMa types for example — though admittedly these are few in number. 
The observed distribution in minimum masses, corrected for completeness for periods up to 
12 years, is shown in Fig. |8l 

The main features seen in this P< 12 yr distribution are (i) a sub-stellar companion mass 
function rising strongly below 10 Mj, (ii) a comparatively flat distribution of stellar compan¬ 
ions, and (iii) a region from ~20-70Mj (corresponding to the brown dwarf regime) where 
relatively few objects are found despite a selection bias in the observations making them 
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easier to detect than pl anets. Note that although objects do exist in this part of the param¬ 


eter space (e.g. 


Jenkins et ah 


20091) , the term brown dwarf des ert was given to 


relative paucity in comparison to planets and stellar objects f lMarcy fc Butler 


righhght the 


200Oj) • Such 


features accord with RV observations elsewhere: the CORALIE, Keck, and Lick surveys all 
report the : same form of sub-s tel lar function while the ‘flat’ stellar distribution mirrors that 
seen in Duqnennov & Ma y or fll99ll. Fig. 11) This similarity comes about despite the fact 


that the various RV surveys work with different samples and operating strategies. The form 
of the distribution of planetary and stellar companions is considered to reflect the differ¬ 
ent formation mechanisms for these two populations (respectively accretion in dissipative 
circumstcllar disks and gravitational instabilities in collapsing cloud fragments) and their 


consequent evolution. The rel atively small number of brown dwar : 


noted elsewhere (for example 


Halbwachs et al. 


2000 


companions has been 


Butler et al. 

2001; 

McCarthy 

2001) 


and may be a re flecti on of a formation mechanism different again from that of stars or plan¬ 
ets, though Armitage & Bonnell ( 2002 ) argue that its existence is a consequence of orbital 
migration of brown dwarfs within an evolving proto-stellar disk. 

The question arises as to what effect a correction for inclination would have on the 
observed distribution. The simplest (crudest) correction is to scale masses up by a factor of 
l/< sinz >. This doesn’t alter the form of the distribution, nor in this case the total number 
of objects that can be regarded as brown dwarf candidates - just those objects out of the 
survey that can be regarded as such. Clearly, a knowledge of the precise form of the mass 
distribution for this (or any) RV survey is precluded until inclination s can be d ete rmined 
accurately via sub-m illi- arcse cond astrometric surveys such as GAIA (ISozzetti et alJ 2001) 


and SIM (jSozzetti et al 


2002]), though it would take a remarkable confluence of inclinations 


for the objects uncovered by the AAPS to alter the underlying distribution of masses. In deed, 


recent attempts have been made to recover the ’true’ mass distribution (Lopez & Jenkins 


2012), with very few of the sub-10 Mj objects moving above this mass limit. 


3.6 Metallicity-Mass Distribution 

One of the most interesting features to emerge from the early study of exoplanets, is_the 


dependence of gas giants to be found orbii 

ting st 

1997; 

Fischer & Valenti 

2005; 

Sousa et al. 

2011) 


accretion scenario for planet formation 


;ars with super-solar metallicities ( Gonzalez. 


20111). This result is a key prediction of the core 


Ida fe Lin 


2004: 


Mordasini et ah 


2012). However, 
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it seems that this bias towards the most meta 


not lower-mass rocky planets ( Udrv et al. 


2007 


•rich stars is on 


Buchhave et al. 


y found for gas giants an d 


2012 


Jenkins et al. 


Therefore, given there is a clear mass dependence as a function of metallicity, it is interesting 
to test what the metallicity distribution looks like for binaries drawn from a representative 
sample. 

In Fig. IH] we show the distribution in the metallicity-mass plane of planets, brown dwarfs, 
and stellar binaries that have been detected mostly by the radial velocity method, with 
a large clutch of the stellar binaries being drawn from the sample of F,G, and K stars 
from Tokovinin et al. (2014). The iron abundances used in this plot were taken from high- 
resolution spectral analysis where possible, generally from the published papers for the 
detected exop lanets, with the Tokovinin et al. primary star metallicities being drawn from 


Casagrande et al. ( 20111 ). The giant planet metallicity bias discussed above is evident here, 


where the sample mean cross-hairs are clearly offset from the sample means at higher masses 
(i.e. above the horizontal dot-dashed line). The brown dwarf and stellar binaries have mean 
values in good agreement with each other, both with sub-solar values, in comparison to the 
exoplanet primary mean distribution that is significantly above the solar value. 



in both of these samples are discussed in each of the works, yet they are large enough and 
have been examined well enough that they can be thought of as good representations of 
their respective held populations. The full Tokovinin et al. sample was not included since we 
wanted to maintain metallicity homogeneity and also we aimed for a direct overlap in orbital 
separation with the exoplanet sample, meaning we only included binary companions with 
orbital periods out to 4 years. For metallicity homogeneity, we cross-matched the sub-4 year 
sample with the ICasagrande et al.1 ( 20111 ) catalogue of metallicities, leaving a total of 874 
binaries or multiple stars, out of the complete 3936 sample. 


As was shown in 


Fischer & Valenti (2005), the distribution of exoplanets follows a power 


law where the percentage of stars with giant planets increases as a function of metallicity. 
After constructing a similar histogram of values to that in Fischer & Valenti, we find a power 
law described by an amplitude of 0.028±0.002 and an index of 2.23±0.09, which we represent 
by the dashed black curve in the plot. Beyond around a solar metallicity, the increase in 






































Spectroscopic Binaries from the AAPS 25 


Companion Mass log(Mf) 


o 

i 


o 

o 


O 

oo 


o 


o 


? p 
x 10 


o 

o 


o 

to 


o 


o 

Oi 


o 

I — 1 
11 1 - 


o 

to 


o 

CO 


o 


• so 

•• • : 

%• o® • 


*.: * 


• 


* X 

*r 


sSu * 


H- 


% % • •• • 




••t" . . *•" 

• ."•Jp-?. ••• • , 

'•*•*:* : 

• • •• • • • 

.;%• . .. 



Figure 9. Metallicity-Mass distribution. Filled circles correspond to exoplanets with 1.5 M .j < M < 5.0 Mj, discovered usin g 


1991 b 


20141) . 


the RV technique, open triangles correspond to binary stars from Halbwachs et all (1201214 and iDuquennov &; Mayor 

filled triangles show our AAPS cand idates, and open circles show the distribution for doubles reported by I Tokovinin __ 

The metallicities were obtained from lCasagrande et alj J201ll) . The horizontal and vertical lines mark the approximate planet- 
brown dwarf boundary and the extreme lower tail of the Tokovinin binary distribution, respectively. The cross-hairs mark the 
sample medians for the planets, AAPS candidates, and the Tokovinin et al. binaries, increasing in mass respectively. 
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[Fe/H] 


Figure 10. The metallicity distribution of exoplanet (upper panel) and binary (lower panel) companions. The best fit model 
distributions are also shown in the plots, where a power-law is used for the exoplanet sample, and a linear function is used to 
model the binary distribution. 


planet hosting fraction rises steeply, possibly accelerating beyond a value of +0.2 dex in 
metallicity flSousa et al.l 12011). 

The distribution of binaries on the other hand is extremely flat across all the metallicity 
rang e, w ithin the uncertaint ies, with the fraction found to be 43+4%, in excellent agreement 


with 


Raghavan et ah ( 2010i ). The best fit weighted linear function is shown in the figure and 


has values for the gradient ( b ) and offset (a) of 0.286+0.011 and -0.099+0.044, respectively. 
These parameters are fairly insensitive to different bin widths, therefore, it is clear that the 
binarity fraction as a function of metallicity is significantly different to the planetary system 
fraction. 


We investigated if there is any dependence on the binary fraction with orbital separation 
by constructing the same distribution on the sample of binaries with orbital periods longer 
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Figure 11. The cumulative fraction of the binary fraction distribution (solid curve) and the exoplanet fraction distribution 
(dashed curve). 


than 4 years, and on the full sample regardless of orbital period. We found no significant 
differences between the distributions, however we do note a small drop in the fraction of 
binaries in the metal-poor regime (-0.8^ [Fe/H] ^-0.5) for the longer period sample. Al¬ 
though it still agrees with a flat distribution within the uncertainties, in the future it may 
be worthwhile to revisit this regime with more binaries in these bins to see if this drop 
in fraction becomes significant, which would indicate there is a dependence of the binary 
fraction with separation as a function of metallicity. 

In Fig. [lT]we show the cumulative fraction between the binary fraction and the exoplanet 
fraction as a function of metallicity. The early rise at low metallicities in the binary fraction 
is apparent, along with the rapid rise at the high metallicity end for the exoplanet fraction. 
We find that the largest disparity between the two populations occurs around 0.15 dex in 
metallicity, a little over solar metallicity, where the planet fraction begins to significantly in¬ 
crease. A Kolmogorov-Smirnov test yields a D-statistic of 0.733 here, relating to a probability 
of the null hypothesis of 2.377xl0~ 2 %. 


Gaoetah 


( 20141 ) show that the fraction of binaries with orbital periods of less than 


1000 days dramatically falls in the metal-rich regime when they split their SDSS and LAM- 
OST FGK star samples up into three different metallicity bins. They find a total binary 
fraction of 43±2%, again in excellent agreement with what we find here, but they calculate 
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binary fractions of 56±5%, 56±3%, and 30±5.7% for metallicities ([Fe/H]) of <-1.1, -1.1 
to <-0.6, and ^-0.6, respectively, indicating a drop in the fraction of binaries in the most 
metal-rich region. Although this is an indication of a higher fraction of binaries near our 
metallicity lower limit, within the uncertainties the distribution is still flat (similar to the 


result reported in 


Raghavan et al 


2010 ). Therefore, we can not confirm if such a binary 


fraction change exists, a change that is also recovered by some models (e.g. 


Machida et al. 


2009 . 


Bate (2014) studied the effects of changes of the metallicity on their star formation 


models, assuming the dust opacity scales linearly with metallicity across a range of metal¬ 
licities from l/3rd solar to three times solar metallicity. He finds no significant changes in 
the multiplicity fraction with metallicity, suggesting gas opacity does not overtly change the 
large-scale properties of star formation from the fragmentation of giant molecular clouds. 
Furthermore, he also finds no dependence with metallicity on the orbital separation of bi¬ 
naries and higher order multiples, as we find here. Given that recent works have shown 
that the cooling times in proto-planetary disks increase as a function of metallicity, meaning 


suppression of disk fragmentation in the super metal-rich regime (e.g. 


Cai et al 


2006), we 


might expect that if FGK star secondaries formed primarily through fragmentation of the 
proto-planetary disk, there would be a strong dependence of the binary fraction with metal¬ 
licity. As this does not appear to be the case, then a flat distribution of binary fraction with 
metallicity suggests that these close binaries predominantly form through fragmentation of 
the giant molecular cloud that also formed the primary star. 


3.7 Summary 

Our target sample of 178 solar-type stars has revealed that ~10 percent are spectroscopic 
binaries. Orbital solutions indicate that two systems potentially have brown dwarf compan¬ 
ions and another two could have eccentricities that place them in the extreme upper tail 
of the eccentricity distribution for binaries with periods less than 1000 days. The systems 
with the largest quantity of data points appear to generate the least robust orbital solutions, 
which could owe to secondary flux contamination of the template spectra, and hence the 
necessity to garner more data to constrain their solutions. When the radial velocity mea¬ 
surement errors are low, yet the Keplerian solutions have low significance, the most likely 
scenario is a multiple-star system. HIPPARCOS astrometry was examined in an attempt 
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to constrain our orbital parameters, however no significant astrometric variation could be 
discerned in the positional residuals. The distribution of companion masses was examined 
for both the binaries and candidate exoplanets detected by the AAPS. For periods up to 
12 years the ‘steep’ planetary and ‘flat’ binary mass distributions mirror those seen by other 
surveys. Over a time scale equivalent to one orbital period of Jupiter, upwards of 30 exo¬ 
planet detections can be expected from our original sample of 178 stars. The discovery of 
these sixteen AAPS binaries from a sample of solar-type stars selected to have no resolvable 
or known SB companions is a reminder that the data for even the relatively bright southern 
stars remain far from complete. 

Finally, analysis of the metallicity-mass plane from planetary companions all the way up 
to stellar companions reveals a stark difference in the mean metallicities of these populations, 
with planets orbiting stars more metal-rich than stellar companions, in general. The fraction 
that host these companions as a function of metallicity is also different, with the binary 
fraction being found to be flat (43±4%) in the metallicity range between -1.0 to 0.6 dex. 
This is in contrast to the fraction of stars hosting giant planetary systems, which has a 
very low fraction until around a solar metallicity, at which point the fraction rises steeply 
following a power law. The flat binary fraction across this wide range of metallicities is in 
agreement with recent hydrodynamical simulations of star formation through fragmentation 
of giant molecular clouds. Such a result suggests this is the dominate formation mechanism 
for FGK-type binaries and not fragmentation of the proto-planetary disk that was left over 
from the formation of the primary star. 
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Table 3. 


HD18907 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s *) 


UCLES 

831.05035 

-616.2 

1.8 



1035.32546 

-245.1 

1.4 



1211.96992 

0.00 

2.3 



1236.91235 

36.6 

2.1 



1383.33260 

203.4 

1.4 



1745.29974 

393.0 

1.6 



1828.11660 

383.8 

2.6 



1856.11391 

393.7 

2.6 



2092.28943 

353.0 

1.6 



2127.24866 

350.0 

2.5 



2710.88687 

-101.4 

2.1 



5101.25580 

-2772.4 

1.3 



5374.32732 

-3044.1 

1.2 

HD25874 

Instrument 

Date (JD - 2450000) 

RV (m s 1 ) 

RV error (m s — *) 


UCLES 

1118.12218 

-1005.7 

1.5 



1473.26127 

-551.9 

1.9 



1526.01361 

-499.2 

1.6 



1630.87612 

-353.4 

1.7 



1830.11789 

-105.3 

2.1 



1920.03798 

0.0 

2.0 



2189.17694 

352.8 

3.0 



2511.23925 

773.4 

4.3 



2594.08090 

882.2 

2.8 



2654.06184 

950.8 

2.2 



2710.89186 

1024.2 

2.4 



6746.88341 

11860.0 

5.2 


HARPS 

5846.82236 

2700.0 

2.1 



5846.82732 

2701.3 

2.1 



5850.82895 

2699.5 

2.1 



5850.83381 

2698.9 

2.1 



5851.81510 

2698.4 

2.1 



5851.82011 

2698.8 

2.1 



5852.81406 

2697.6 

2.1 



5852.81897 

2697.2 

2.1 



5984.50942 

2669.3 

2.0 



5984.51441 

2669.2 

2.1 



5986.53569 

2668.6 

2.0 



5986.54065 

2668.6 

2.0 



6298.59813 

2585.0 

2.1 



6298.60323 

2585.8 

2.1 



6301.61764 

2587.3 

2.1 



6301.62250 

2585.5 

2.1 



6746.88000 

2453.6 

5.2 

HD26491 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s — *) 


UCLES 

831.07678 

-2115.6 

1.5 



1118.13648 

-784.2 

1.6 



1473.26423 

0.0 

2.0 



1525.97709 

70.9 

1.8 



1527.06295 

65.8 

1.8 



1745.32814 

274.1 

1.9 



1920.04374 

349.3 

3.3 



2127.30280 

439.6 

3.0 



2710.89727 

413.5 

1.8 



5101.20238 

-697.0 

1.3 



5845.24838 

-1236.9 

2.9 



5899.07115 

-1278.6 

2.3 
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Table 3. (cont.) 


HD39213 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s x ) 


UCLES 

1118.21233 

13.2 

5.9 



1121.17619 

12.6 

7.6 



1157.15034 

-122.2 

8.1 



1212.99066 

-361.0 

8.0 



1214.09191 

-381.4 

7.5 



1274.86986 

-553.7 

7.5 



1275.87068 

-545.8 

8.2 



1387.33431 

-809.0 

10.9 



1683.85561 

772.0 

10.6 



1828.19855 

1638.4 

9.1 



1830.04517 

1657.9 

12.6 



1856.2322 

1681.6 

9.5 



2127.33354 

1238.2 

19.1 



2751.87683 

-777.0 

6.8 

HD42024 

Instrument 

Date (JD - 2450000) 

RV (m s" 1 ) 

RV error (m s -1 ) 


UCLES 

1118.2294 

-2859.2 

8.0 



1157.16631 

3573.5 

2.9 



1213.02311 

-1299.5 

2.8 



1214.10000 

-1052.7 

2.7 



1275.87845 

-2848.1 

3.5 



1411.31272 

-1804.3 

5.8 



1414.30077 

-2188.3 

2.5 



1473.23953 

1382.8 

2.8 



1530.15077 

2200.9 

3.1 



1630.95238 

0.0 

3.2 



1920.06937 

3570.0 

3.3 



1983.93173 

964.8 

3.4 



2060.84424 

1201.9 

2.9 

HD64184 

Instrument 

Date (JD - 2450000) 

RV (m s 1 ) 

RV error (m s -1 ) 


UCLES 

946.85278 

0.0 

2.7 



1236.06303 

-8217.3 

2.8 



1630.99829 

2005.9 

2.0 



1683.88873 

-1894.8 

1.9 



1920.11122 

13652.9 

2.8 



1983.99441 

-5721.9 

3.5 



2009.04858 

12902.0 

2.9 



2751.98869 

-5134.5 

2.4 
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Table 3. (cont.) 


HD120690 Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s 1 ) 

UCLES 

917.11028 

81.0 

3.7 


1236.24332 

-4831.4 

3.1 


1275.12458 

-4967.8 

4.0 


1382.94887 

-4855.0 

2.5 


1630.26344 

-2914.8 

2.7 


1683.06432 

-2373.5 

2.9 


1717.96127 

-2028.6 

2.8 


1743.00961 

-1771.8 

2.8 


1920.24338 

0.0 

3.3 


1984.14882 

580.2 

4.2 


2060.96853 

1224.4 

3.4 


2092.97458 

1483.7 

3.1 


2748.21844 

5322.1 

5.5 


6793.06682 

6289.8 

3.2 

HARPS 

5983.85998 

-897.3 

2.0 


5983.86440 

-898.1 

2.1 


5984.84268 

-896.4 

2.1 


5984.84727 

-893.5 

2.1 


5985.83524 

-889.2 

2.1 


5985.83971 

-889.3 

2.1 


5986.85772 

-877.0 

2.1 


5986.86206 

-876.2 

2.1 


6042.61245 

-472.6 

2.1 


6042.61687 

-473.0 

2.1 


6042.76914 

-472.1 

2.1 


6042.77348 

-473.0 

2.1 


6046.73467 

-446.5 

2.1 


6046.73870 

-443.9 

2.1 


6047.75428 

-439.1 

2.1 


6047.75862 

-438.9 

2.1 


6048.78013 

-435.9 

2.1 


6048.78467 

-434.6 

2.1 


6300.85104 

1104.6 

2.1 


6300.85551 

1105.5 

2.1 


6375.78422 

1515.5 

2.0 


6375.78865 

1516.6 

2.0 


6376.80760 

1515.1 

2.1 


6376.81202 

1515.6 

2.1 
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Table 3. (cont.) 


HD121384 Instrument Date (JD - 2450000) RV (m s 1 ) RV error (m s -1 ) 


UCLES 

831.23238 

-761.1 

1.2 


915.11284 

1173.1 

2.2 


1002.94245 

-119.4 

2.3 


1212.28329 

-9052.3 

3.6 


1236.25493 

289.7 

2.8 


1383.96012 

-9086.6 

2.1 


1385.85181 

-14785.4 

2.3 


1386.86423 

-19043.8 

2.2 


1410.87119 

0.0 

3.9 


1413.85641 

211.3 

3.5 


1631.26748 

1160.1 

2.4 


1682.99133 

992.2 

3.0 


1684.06855 

980.7 

2.8 


1706.07667 

451.6 

3.3 


1717.83328 

-155.9 

3.1 


1742.88968 

-13609.7 

2.8 


1743.88168 

-17735.1 

2.7 


1919.25971 

-7795.4 

5.3 


1984.16157 

1114.1 

2.3 


2010.26683 

1205.2 

1.8 


2061.01149 

535.1 

2.7 


2091.91315 

-2958.7 

3.6 


2127.89869 

145.3 

2.4 


2129.92347 

270.0 

2.3 


2752.11488 

1007.9 

1.1 


5669.18415 

-5277.9 

0.8 

HD 131923 Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error(m s _1 ) 

UCLES 

831.26615 

8714.9 

1.8 


915.14155 

8120.8 

1.9 


1236.25951 

6268.9 

2.5 


1383.98096 

5644.7 

2.8 


1683.01896 

4640.8 

2.3 


1718.05774 

4541.6 

3.2 


1984.19127 

3894.3 

2.2 


2187.86419 

3481.1 

2.8 


2711.31697 

2566.7 

3.5 


2748.22162 

2505.1 

3.0 


5013.85913 

20.8 

1.2 


5017.87936 

18.0 

1.1 


5018.99094 

13.0 

1.0 


5019.96696 

14.7 

1.1 


5020.90813 

6.1 

1.4 


5021.91260 

7.3 

1.4 


5023.90683 

-0.0 

1.5 


5029.89593 

-3.4 

1.2 


5030.85082 

-11.1 

1.4 


5031.96501 

-2.9 

1.1 


5032.99508 

-1.5 

1.3 


5036.89032 

-13.7 

1.3 


5037.88830 

-17.3 

1.6 


5043.92813 

-24.9 

1.3 


5044.87581 

-29.0 

1.3 


5045.89427 

-25.5 

1.1 


5048.90380 

-23.05 

1.2 


5049.88025 

-26.0 

1.7 


6793.09817 

-1791.3 

2.5 
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Table 3. (cont.) 


HD131923 (cont.) 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error(m s *) 


HARPS 

5983.91673 

-1033.9 

2.1 



5983.91976 

-1033.0 

2.1 



5983.92322 

-1036.3 

2.0 



5984.89657 

-1022.5 

2.1 



5984.89980 

-1023.3 

2.1 



5984.90305 

-1023.6 

2.1 



5985.89276 

-1023.4 

2.0 



5985.89603 

-1026.9 

2.0 



5985.89926 

-1024.6 

2.0 



6042.69551 

-1005.8 

2.1 



6042.69889 

-1006.0 

2.1 



6042.70204 

-1007.8 

2.1 



6042.70610 

-1007.7 

2.1 



6042.70944 

-1004.9 

2.1 



6042.71259 

-1008.4 

2.1 



6046.77564 

-1006.8 

2.1 



6046.77888 

-1005.3 

2.1 



6046.78219 

-1006.6 

2.1 



6047.80434 

-1005.1 

2.1 



6047.80758 

-1006.9 

2.1 



6047.81077 

-1005.2 

2.1 



6048.82939 

-1011.2 

2.1 



6048.83265 

-1010.8 

2.1 



6048.83586 

-1009.0 

2.1 



6165.48087 

-961.2 

2.0 



6165.48409 

-960.8 

2.0 



6165.48724 

-960.9 

2.0 



6299.84138 

-939.4 

2.1 



6299.84456 

-940.9 

2.1 



6300.87569 

-930.3 

2.1 



6300.87892 

-931.1 

2.1 



6300.88222 

-929.2 

2.1 



6301.86354 

-927.0 

2.1 



6301.86687 

-925.0 

2.1 



6301.87007 

-926.2 

2.1 



6375.83888 

-927.1 

2.0 



6375.84214 

-926.0 

2.0 



6375.84538 

-927.2 

2.0 



6376.86779 

-928.7 

2.1 



6376.87102 

-927.8 

2.1 



6376.87431 

-928.3 

2.1 



6377.79462 

-932.6 

2.1 



6377.79798 

-930.0 

2.1 



6377.80119 

-930.7 

2.1 

HD145825 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error(m s — *) 


UCLES 

915.18178 

-348.0 

2.2 



1002.04574 

-448.0 

3.0 



1382.97383 

-534.0 

1.9 



1630.27976 

-375.7 

2.0 



1683.04691 

-330.1 

2.1 



1718.09632 

-310.8 

2.2 



1742.99728 

-287.7 

2.1 



1766.89725 

-269.1 

1.8 



1984.22396 

-79.4 

2.3 



2060.98194 

0.0 

2.0 



2091.94491 

35.1 

2.1 



2126.92735 

56.7 

2.4 



2711.31479 

544.4 

3.0 



2748.21470 

564.8 

2.7 



5669.19403 

1646.1 

1.8 



6088.20721 

1610.5 

2.6 



6765.20219 

1103.0 

1.2 
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Table 3. (cont.) 


HD145825 (cont.) 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s *) 


HARPS 

6042.71940 

277.4 

2.1 



6042.72434 

278.2 

2.1 



6046.82852 

271.8 

2.1 



6046.83381 

271.3 

2.1 



6047.87475 

275.3 

2.1 



6047.87965 

275.6 

2.1 



6048.88054 

276.2 

2.1 



6048.88548 

276.7 

2.1 



6164.53742 

208.8 

2.1 



6164.54216 

208.7 

2.1 



6165.51908 

207.6 

2.1 



6165.52394 

206.5 

2.1 



6375.90768 

36.4 

2.0 



6375.91245 

36.6 

2.0 



6376.89108 

35.6 

2.1 



6376.89585 

36.1 

2.1 



6377.86316 

34.3 

2.1 



6377.86788 

34.4 

2.1 



6557.50361 

-128.5 

2.4 



6557.50834 

-130.6 

2.4 



6558.48109 

-131.0 

2.1 



6558.48582 

-130.1 

2.1 



6559.49101 

-128.4 

2.1 



6559.49573 

-129.1 

2.1 



6560.48828 

-126.6 

2.1 



6560.49300 

-127.6 

2.1 

HD150248 

Instrument 

Date (JD - 2450000) 

RV (m s ^) 

RV error (m s _1 ) 


UCLES 

917.26297 

-506.1 

3.0 



1384.03661 

-413.2 

2.2 



1683.07465 

-229.2 

2.3 



1684.12307 

-232.1 

2.3 



1743.03057 

-164.5 

2.1 



1766.91499 

-142.2 

2.2 



1767.94821 

-134.0 

2.4 



1827.89080 

-62.5 

2.5 



1828.87861 

-54.9 

2.2 



1984.23953 

259.2 

2.7 



2061.09080 

541.6 

2.3 



2091.96333 

670.8 

2.0 



2126.95164 

865.5 

2.0 



2389.21141 

3448.5 

1.0 



2711.31186 

612.6 

2.2 



2748.22495 

499.7 

3.1 


HARPS 

6042.74497 

244.9 

2.0 



6046.86431 

231.2 

2.0 



6047.90016 

228.8 

2.0 



6048.89360 

226.6 

2.0 



6164.57497 

-24.5 

2.1 



6165.55598 

-26.3 

2.0 
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Table 3. (cont.) 


HD156274B Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s 1 ) 

UCLES 

915.25103 

-100.2 

1.8 


1236.29045 

-93.6 

1.9 


1631.31256 

-75.6 

1.4 


1683.08190 

-77.1 

1.9 


1684.13026 

-74.9 

1.9 


1718.10916 

-69.6 

1.6 


1743.97223 

-69.9 

1.5 


1766.92942 

-72.9 

1.7 


1767.96359 

-73.8 

1.6 


1984.25459 

-66.3 

1.8 


2061.10543 

-58.2 

1.7 


2093.06873 

-64.4 

1.7 


2127.96561 

-68.4 

1.9 


2151.96653 

-62.5 

1.8 


2186.90647 

-57.3 

1.4 


2189.90668 

-58.9 

1.3 


2359.21366 

-51.1 

1.4 


2387.16082 

-54.1 

1.4 


2388.20299 

-55.3 

1.5 


2422.14228 

-48.9 

1.3 


2455.03140 

-50.8 

1.6 


2509.95469 

-49.0 

1.7 


2510.92182 

-54.4 

1.5 


2745.25045 

-47.2 

1.6 


2746.30657 

-46.1 

1.4 


2752.16529 

-45.4 

1.6 


2784.11807 

-48.4 

1.5 


2858.94418 

-47.6 

1.6 


2942.92355 

-40.3 

1.6 


3041.27889 

-38.1 

1.7 


3042.29099 

-45.0 

1.7 


3214.97664 

-26.2 

1.5 


3245.00697 

-31.9 

1.7 


3280.90584 

-36.0 

1.5 


3483.23536 

-25.1 

0.9 


3486.15507 

-24.5 

0.7 


3507.14609 

-26.5 

0.8 


3516.11359 

-32.2 

1.0 


3517.15245 

-28.0 

0.9 


3520.20139 

-22.4 

0.9 


3521.18310 

-24.9 

0.8 


3523.11361 

-21.8 

0.8 


3569.02820 

-23.7 

0.8 


3627.88388 

-32.6 

1.1 


3628.91553 

-25.2 

0.8 


3631.88150 

-25.8 

0.7 


3842.22656 

-20.6 

0.9 


3938.05195 

-12.5 

0.8 


4010.91205 

-20.1 

0.7 


4226.20607 

-7.9 

1.0 


4255.00536 

-8.0 

0.9 


4336.01844 

-9.2 

1.1 


4550.26395 

-3.0 

1.8 


5013.93467 

1.4 

0.9 


5014.99159 

1.8 

2.0 


5017.91353 

2.3 

1.0 


5019.04922 

2.4 

1.0 


5020.07329 

5.4 

1.0 


5020.86256 

-0.2 

1.0 


5021.99621 

4.5 

0.9 
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Table 3. (cont.) 


HD156274B (cont.) Instrument Date (JD - 2450000) RV (m s x ) RV error (m s x ) 


UCLES 

5023.00447 

1.6 

1.1 


5024.00843 

7.2 

0.9 


5027.06874 

6.2 

0.9 


5029.96115 

3.1 

1.0 


5030.12974 

5.3 

0.9 


5030.91509 

1.8 

0.8 


5032.04541 

4.7 

0.9 


5032.96709 

5.4 

0.9 


5036.96946 

9.2 

0.8 


5039.01974 

13.0 

0.8 


5040.00608 

10.5 

1.4 


5043.97142 

8.1 

0.7 


5044.96994 

5.0 

0.9 


5046.86756 

5.4 

0.6 


5047.97102 

3.8 

0.9 


5049.06488 

1.7 

0.7 


5049.94976 

1.6 

0.9 


5050.94561 

-3.0 

0.9 


5051.91615 

-1.7 

0.9 


5052.94001 

-2.1 

1.6 


5053.94417 

2.9 

1.3 


5054.94608 

3.2 

0.9 


5056.03535 

4.0 

0.8 


5056.99612 

6.4 

0.8 


5057.99625 

4.9 

0.8 


5058.98884 

6.5 

0.8 


5076.03618 

-0.7 

1.1 


5100.89570 

-1.3 

1.6 


5110.92624 

2.0 

1.1 


5310.15913 

13.1 

0.9 


5374.08325 

16.5 

2.0 


5403.07787 

16.4 

1.1 


5996.20104 

34.9 

1.1 


6088.21396 

40.8 

1.7 


6767.15356 

52.1 

1.0 

HARPS 

3039.89653 

-124.3 

2.1 


3039.89901 

-125.4 

2.1 


3039.90142 

-131.1 

2.1 


3306.48142 

-118.9 

2.0 


3306.48346 

-120.3 

2.0 


3306.48557 

-119.5 

2.0 


3429.78850 

-115.1 

2.0 


3429.79210 

-115.3 

2.0 


3429.79577 

-116.8 

2.0 


3429.85697 

-116.4 

2.0 


3429.86064 

-116.5 

2.0 


3429.86443 

-116.5 

2.0 


3521.77104 

-114.7 

2.0 


3521.77480 

-116.0 

2.0 


3521.77880 

-115.5 

2.0 


3523.11000 

-107.5 

0.8 


3537.78318 

-115.9 

2.0 


3537.78700 

-114.8 

2.1 


3537.79085 

-115.5 

2.1 


3600.47566 

-111.7 

2.0 


3600.47931 

-112.5 

2.0 


3600.48326 

-112.0 

2.0 


3808.76299 

-108.9 

2.0 
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Table 3. (cont.) 


HD156274B (cont.) 

Instrument 

Date (JD - 2450000) 

RV (m s” 1 ) 

RV error (m s *) 


HARPS 

3808.76486 

-109.5 

2.0 



3808.76669 

-109.3 

2.0 



3808.78437 

-109.1 

2.1 



3808.78615 

-108.3 

2.1 



3808.78793 

-108.9 

2.0 



3808.80520 

-109.3 

2.0 



3808.80703 

-109.1 

2.0 



3808.80886 

-107.5 

2.0 



3828.81348 

-107.0 

2.0 



3828.81529 

-106.9 

2.0 



3828.81703 

-107.1 

2.0 



3873.65716 

-107.2 

2.0 



3873.65894 

-108.3 

2.0 



3873.66069 

-107.3 

2.0 



3877.67621 

-108.4 

2.0 



3877.67807 

-108.8 

2.0 



3877.67983 

-107.9 

2.0 



3896.75731 

-103.2 

2.1 



3896.75908 

-103.0 

2.1 



3896.76077 

-103.7 

2.1 



3915.84989 

-106.3 

2.0 



3915.85169 

-107.8 

2.0 



3915.85343 

-105.5 

2.0 



3988.62757 

-104.6 

2.1 



3988.62934 

-103.3 

2.1 



3988.63117 

-104.7 

2.0 



4161.80181 

-100.4 

2.1 



4161.80356 

-99.6 

2.1 



4161.80531 

-98.8 

2.1 



4192.86533 

-100.3 

2.0 



4192.86708 

-100.0 

2.0 



4192.86883 

-99.5 

2.0 



4203.89316 

-100.0 

2.0 



4203.89489 

-102.0 

2.0 



4203.89665 

-101.3 

2.0 



4246.68353 

-96.1 

2.1 



4246.68528 

-96.7 

2.1 



4246.68703 

-96.7 

2.1 



4290.69783 

-97.1 

2.2 



4290.69997 

-96.3 

2.2 



4290.70167 

-96.7 

2.3 



4312.59589 

-95.0 

2.0 



4312.59765 

-96.1 

2.1 



4312.59944 

-95.8 

2.0 



4339.53113 

-94.1 

2.0 



4339.53288 

-95.0 

2.0 



4339.53466 

-94.9 

2.0 



4350.55236 

-96.1 

2.0 



4350.55412 

-94.4 

2.1 



4350.55587 

-95.8 

2.1 

HD158783 

Instrument 

Date (JD - 2450000) 

RV (m s 1 ) 

RV error (m s — *) 


UCLES 

915.25887 

-1160.7 

5.3 



1384.06035 

-891.6 

1.8 



1386.87073 

-871.0 

1.8 



1410.89183 

-828.4 

2.2 



1413.89204 

-821.8 

1.7 



1630.30834 

-320.9 

2.0 



1683.08959 

-182.0 

2.2 



1684.13587 

-185.8 

2.2 
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Table 3. (cont.) 


HD158783 (cont.) 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s *) 


UCLES 

1706.10237 

-107.5 

2.1 



1718.11425 

-82.0 

2.0 



1742.90540 

-6.2 

2.0 



1743.91951 

-9.8 

2.4 



1745.05091 

0.8 

2.1 



1827.90831 

251.1 

2.2 



1828.88293 

246.5 

2.5 



1829.88574 

251.8 

2.8 



1984.25818 

725.0 

3.0 



2009.18580 

826.7 

2.7 



2061.10905 

985.6 

2.3 



2091.97455 

1090.1 

2.1 



2126.96948 

1188.2 

1.9 



2711.30821 

2693.2 

1.4 



6793.23833 

1587.9 

1.2 



6794.20468 

1598.6 

1.0 

HD162255 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s _1 ) 


UCLES 

1002.05244 

6218.0 

4.7 



1382.99601 

1036.8 

2.0 



1410.96927 

-11071.1 

2.9 



1630.29429 

13912.1 

2.1 



1684.17416 

16122.8 

2.6 



1718.13792 

0.0 

1.9 



1742.94071 

-3974.7 

2.2 



1744.02728 

-6307.1 

1.9 



1766.96237 

1667.6 

2.1 



1827.92694 

16187.1 

2.4 



1984.28228 

-7283.5 

2.9 



2061.13815 

12969.8 

2.0 



2091.99501 

-14472.3 

2.3 



2127.03219 

-4974.4 

1.6 

HD169586 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s _1 ) 


UCLES 

917.30325 

-6740.6 

6.5 



1120.89231 

-5470.6 

4.6 



1683.11045 

-892.4 

5.0 



1744.01838 

-152.0 

5.3 



1745.08961 

-160.5 

5.0 



1766.96649 

0.0 

5.4 



2061.14242 

3385.1 

7.1 



2711.31987 

863.5 

11.2 



2748.20956 

-508.7 

10.2 



5374.13437 

6244.4 

7.2 

HD1T5345 

Instrument 

Date (JD - 2450000) 

RV (m s- 1 ) 

RV error (m s _1 ) 


UCLES 

1119.92272 

-57.5 

3.0 



1387.19787 

3319.5 

3.0 



1410.95038 

1550.7 

3.5 



1472.91236 

-3628.9 

3.2 



1683.12996 

4575.2 

4.3 



1706.15776 

2859.0 

8.9 



1718.14493 

1940.2 

3.3 



1742.95621 

73.7 

4.1 



1743.99665 

0.0 

3.8 



1745.09575 

-78.5 

4.4 



1827.94994 

-9067.6 

3.5 



1984.29080 

5441.5 

5.1 



2061.15604 

-367.2 

5.0 



2092.01501 

-3053.1 

4.1 



2127.06057 

-7052.3 

3.1 



2748.20573 

-6569.9 

5.3 






